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We explore the possibility of using the -B-type polarization of the Cosmic Microwave Background 
to map the large-scale structures of the Universe taking advantage of the lens effects on the CMB 
polarization. The functional relation between the B component with the primordial CMB polariza- 
tion and the line-of-sight mass distribution is explicited. Noting that a sizeable fraction (at least 
40%) of the dark halo population which is responsible of this effect can also be detected in galaxy 
weak lensing survey, we present statistical quantities that should exhibit a strong sensitivity to this 
overlapping. We stress that it would be a sound test of the gravitational instability picture, inde- 
pendent on many systematic effects that may hamper lensing detection in CMB or galaxy survey 
alone. Moreover we estimate the intrinsic cosmic variance of the amplitude of this effect to be less 
than 8% for a 100 deg 2 survey with a 10' CMB beam. Its measurement would then provide us with 
an original mean for constraining the cosmological parameters, more particularly, as it turns out, 
the cosmological constant A. 
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I. INTRODUCTION 

In the new era of precision cosmology we are enter- 
ing in, the forthcoming experiments will provides us 
with accurate data on Cosmic Microwave Background 
anisotropics This should lead to accurate determina- 
tions of the cosmological parameters, provided the large- 
scale structures of the Universe indeed formed from grav- 
itational instabilities of initial adiabatic scalar pertur- 
bations. It has been soon realized however that even 
with the most precise experiments, the cosmological pa- 
rameter space is degenerate when the primary CMB 
anisotropies alone are considered Complementary 
data, that may be subject to more uncontrollable sys- 
tematics are thus required, such as supernovae surveys 
H (but see 0) or constraints derived from the large- 
scale structure properties. Among the latter, weak lens- 
ing surveys are probably the safer ||, but still have not 
yet proved to be accurate enough with the present day 
observations. 

Secondary CMB anisotropies (i.e. induced by a subse- 
quent interaction of the photons with the mass or matter 
fluctuations) offer opportunities for raising this degener- 
acy. Lens effects |(J are particularly attractive since they 
are expected to be one of the most important. They also 
are entirely driven by the properties of the dark matter 
fluctuations, the physics of which involve only gravita- 
tional dynamics, and are therefore totally controlled by 
the cosmological parameters and not by details on galaxy 
or star formation rates. More importantly an unambigu- 
ous detection of the lens effects on CMB maps would be a 
precious confirmation of the gravitational instability pic- 
ture. Methods to detect the lens effects on CMB maps 
have been proposed recently. High order correlation func- 



tions Qj], peak ellipticities || or large scale lens induced 
correlators ||] have been proposed for detecting such ef- 
fects. All of them are however very sensitive to cosmic 
variance since lens effect is only a sub-dominant alter- 
ation of the CMB temperature patterns. The situation 
is different when one considers the polarization proper- 
ties. The reason is that in standard cosmological models 
temperature fluctuations at small scale are dominated by 
scalar perturbations. Therefore the pseudo-scalar part, 
the so called B component, of the polarization is negligi- 
ble compared to its scalar part (the E component) and 
can only be significant when CMB lens couplings are 
present. This mechanism has been recognized in earlier 
papers []l0| , pl]| . The aim of this paper is to study sys- 
tematically the properties of the lens induced B field and 
uncover its properties. 

In section ||, we perturbatively compute the lens effect 
on the CMB polarization E and B field. This first order 
equation is illustrated by numerical experiments. Possi- 
bility of direct reconstruction of the projected mass dis- 
tribution is also examined. As it has already been noted 
a significant fraction of the potential wells that deflect 
the CMB photons can actually be mapped in local weak 
lensing surveys |^2| , H - This feature has been considered 
so far in relation to the CMB temperature fluctuations. 
We extend in Section QI these studies to the CMB po- 
larization exploiting the specificities of the field found in 
previous section. In particular we propose two quanti- 
ties that can be built from weak lensing and Cosmic 
Microwave Background polarization surveys, the average 
value of which does not vanish in presence of CMB lens 
effects. Compared to direct analysis of the CMB polar- 
ization, such tools have the joint advantage of being less 
sensitive to systematics -systematic errors coming from 
CMB mapping on one side and weak lensing measure- 
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ment on the other side have no reason to correlate!- and 
so emerge even in presence of noisy data, and of being an 
efficient probe of the cosmological constant. Indeed the 
expected amplitude of correlation is directly sensitive to 
the relative length of the optical bench, from the galaxy 
source plane to the CMB plane, which is mainly sensitive 
to the cosmological constant. Filtering effects and cosmic 
variance estimation of such quantities are considered in 
this section as well. 



II. LENS EFFECTS ON CMB POLARIZATION 

A. First order effect 

Photons emerging from the last scattering surface are 
deflected by the large scale structures of the Universe 
that are present on the line-of-sights. Therefore photons 
observed from apparent direction a must have left the 
last scattering surface from a slightly different direction, 
a + where £ is the lens induced apparent displace- 
ment at that distance. The displacement field is related 
to the angular gradient of the projected gravitational po- 
tential. In the following, the lens effect will be described 
by the deformation effects it induces, encoded in the am- 
plification matrix, 



so observed quantity will not take into account any other 
secondary effect. It is very important at this point to 
note that the lensing effect does not produce any polar- 
ization nor rotate the Stokes parameter. In this regime 
its effect reduces to a simple deformation of the polar- 
ization patterns, similar to the temperature maps. This 
is the mechanism by which the geometrical properties of 
the polarization field are changed. 

To see that we have to consider the electric (E) and 
magnetic (B) components instead of the Stokes param- 
eters. At small angular scales (we assume that a small 
fraction of the sky can be described by a plane), these 
two quantities are defined as, 
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The lens effect affects the local polarization just by mov- 
ing the apparent direction of the line of sight jui). Thus, 
if we use the Stokes parameters Q and U to describe the 
local polarization vector, 



P = 



U 



we can relate the observed polarization P to the primor- 
dial one by the relation 



Q(a) = Q(a + i), U(a) = U{a + £). 



(3) 



From now on we will denote x an observed quantity and 
x the primordial one. ~d' = a + £ is the sky coordi- 
nate system for the observer, therefore the amplification 
matrix A is also the Jacobian of the transformation be- 
tween the source plane and the image plane. We will 
restrain here our computation to the weak lensing effect 



This fields reflect non-local geometrical properties of the 
polarization field. The electric component accounts for 
the scalar part of the polarization and the magnetic one, 
the pseudo-scalar part: by parity change E is conserved, 
whereas B sign is changed. As it has been pointed out in 
previous papers ||l0| , |ll| , |l4| , lens effects partly redistribute 
polarization power in these two fields. 

We explicit this latter effect in the weak lensing regime 
where distortions, k and 7 j components are small. This 
is indeed expected to be the case when lens effects by 
the large-scale structures are considered, for which the 
typical value of the convergence field n is expected to 
be ~ 1% at 1 degree scale. The leading order effect is 
obtained by simply pluging (||) in ([|) and by expanding 
the result at leading order in £ , k, and 7 . Noting that 
(these calculations are very similar to those done in 
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we can write a perturbation description of the lensing 
effect on electric and magnetic components of the polar- 
ization. At leading order one obtains: 

AE = AE + CdiAE - 2nAE 

-2 < 5. iJ ( 7 l AP^+ 7 ; fc P^) + 0( 7 2 ) 

AB = AB + Cd t AB - 2nAB (6) 
-2e^(f AP^+ 7 ; fc P^)+0( 7 2 ), 

Where we used the fact that AX = AX + <f <%AA at 
the leading order. The formulas for E and B are alike. 
The only difference stands in the 5ij and Uj (the lat- 
ter is the totally antisymmetric tensor, en = e 2 2 = 0, 
612 = —£21 = 1) that reflects the geometrical properties 
of the two fields. The first three terms of each of these 
equations represent the naive effect: the lens induced de- 
formation of the E or B fields. This effect is comple- 
mented by an enhancement effect (respectively k AE and 



5«r. 






d l X 










didjX 


= d^X-(6>*+&(5 l j +Q 






(2) 







2 



k AB) and by shear-polarization mixing terms. The latter 
effects consist in two parts. One which we will call the A- 
term that couples the shear with second derivative of the 
polarization field. The other one, hereafter the V-term, 
mixes gradient of the shear and polarization. Although 
terms like V have been neglected in similar computations 
[jl3| we cannot do that here a priori. We will indeed show 
later that these two terms have similar amplitudes. 

One consequence of standard inflationary models on 
CMB anisotropies is the unbalanced distribution of power 
between the electric (E) and magnetic (B) component 
of its polarization. Adiabatic scalar fluctuations do not 
induce B-type polarization and they dominate at small 
scales over the tensor perturbations (namely the gravity 
waves). So, even though gravity waves induce E and B 
type polarization in a similar amount, primary CMB 
sky is expected to be completely dominated by E type 
polarization at small scales. Then for this class of models 
the actual magnetic component of the polarization field 
is generated by the corrective part of eq. (||) , 
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This result extends the direct lens effects described in 
Benabed & Bernardeau filll who focused their analysis 
on the lens effect due to the discontinuity of the polar- 
ization field in case of cosmic strings. Previous studies of 
the weak lensing effect on CMB showed that with lens- 
ing, the B component becomes important at small scales 
[jTJI . We obtain here the same result but with a different 
method; eq. (Q) means that the polarization signal P is 
redistributed by the lensing effect in a way that breaks 
the geometrical properties of the primordial field. Note 
here that it is mathematically possible to build a shear 
field that preserves these geometrical properties and that 
does not create any B signa l at s mall scales. We will dis- 
cuss this problem in Sec. QC. It also means that B 



directly reflects the properties of the shear map. We will 
take advantage of this feature to probe the correlation 
properties of B with the projected mass distribution in 
next sections. 



B. Lens-induced B maps 

We show examples of lens induced B maps. These 
maps have been calculated using "CMBSlow" code de- 
veloped by A. Riazuelo (see JT7J]) to compute primordial 
polarization maps (we use realizations of standard CDM 
model to illustrate lens effects) . Then various shear maps 
are applied. We present both true distortions, (obtained 
by Delaunay triangulation used to shear the Q and U 
fields), and the first order calculations given by eq. ([?]). 

Fig. [l] presents the shear effect induced by an isother- 
mal sphere with finite core radius (and the lens edges 
have been suppressed by an exponential cutoff to min- 
imize numerical noise). The agreement between true 




FIG. 3. The C e of a 100 square degree B map. The solid 
line is the full first order approximation formula. The dotted 
line gives the contribution of the V-term. The dash-dotted 
one represents the A-term. The latter is dominant at small 
is, around I = 1000, that is to say for structures around 10'. 
The V contribution gives birth to smaller structures in the 
1 ~ 2 arc-minute range. 



distortion (central panel ) and first order formula (right 
panel) is good. However, a close examination of the maps 
reveals that some structures in the true map are slightly 
wider than their counterparts in the first order map. This 
error is more severe in the center, where the distortion 
is bigger, which is to be expected since the limits of the 
validity region of first order calculations are reached. 

Fig. H shows the B field induced by a realistic dis- 
tortion. We use second order Lagrangian dynamics |2l]] 
to create a 2.5 x 2.5 degree map that mimics a realistic 
projected mass density up to z = 1000 and used its grav- 
itational distortion to compute a typical weak lensing- 
induced B map. Again we compare the exact effect (i.e. 
left panel where Delaunay triangulation is used) and the 
first order formula (middle panel). Right panel shows 
the difference between the two maps. It reveals the loca- 
tions where the two significantly disagree. In fact most 
disagreements are due to slight mismatch of the B patch 
positions, which lead to dipole like effects in this map. 

We also show here a comparison of the two parts of 
the first order formula eq. (^) in order to see which of 
the A or V terms dominates. It would be more com- 
fortable if one of the two terms was dominant, however, 
Fig. H shows that it is not the case. Even if the A- 
term dominates at low (< 1000) i, it is only twice bigger 
than V-one at this scale. The inverse is true for higher 
(3000 ~ 5000) is. This can be seen by looking at Fig. | 
where we show the relative amplitudes of the A and V 
contributions. The A part gives birth to large patches 
(around 10') while V panel shows a lot more of small 
features. 
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FIG. 1. Lens effect induced by a large isothermal sphere with finite core radius. The k map of the lens is shown on left panel. 
The primordial E sky is presented in the middle left panel. It has been generated for a Slo = 0.3, A = 0.7 model, without tensor 
modes. The middle right panel displays the true reconstructed AB field in a 4.5 x 4.5 degree map and the right panel shows 
the first order approximation. Note that the rosette-like shape the eye seams to catch in B fields is a numerical coincidence 
and has no special significance. 




FIG. 2. The effect of a realistic weak lensing field on B. 2.2 x 2.2 degree survey with 1.8' resolution. The left panel shows 
exact distortion obtained by Delaunay triangulation. The middle one, the first order formula result, and the right gives the 
difference between the two. The three panels share the same color table. The mean amplitude in the difference map is about 
3 times smaller. 



4 



C. Direct reconstruction — Kernel problem 

The fact that the observable B is at leading order pro- 
portional to the weak lensing signal invites us to try a 
direct reconstruction, similar to the lensing mass recon- 
struction. In fact, we can write 



AB = -2e 



k\ _ 



= F[ 7 ] 



(8) 



and our reconstruction problem becomes an inversion 
problem for the operator F. Unfortunately, one can prove 
that this problem has no unique solution. It is due to the 
fact that F admits a huge kernel, in the sense that, given 
a polarization map, there is a wide class of shear fields 
that will conserve a null B polarization. The demonstra- 
tion of this property is sketched in the following. 

Since the unlensed polarization is only electric in our 
approximation, we can describe it by the Laplacian of a 
scalar field ; 



E = Aip so 



Q = {dl - dl) <p 



(9) 



U = 2 d x d y (f 

The same holds for the shear and convergence fields 

« = 7i = \ {dl ~ dl) 4 , 72 = d x 8 y xl>. (10) 

Thus we need to know, for a given (p field, whether there 
is any ip that fulfills the equation 

72 AQ - 7i AC/ + 9 l72 d l Q - 8 i7l d l U = 0. (11) 

(p and ip can be written as polynomial decompositions 

<fi(x,y) = ^ a n ix n y l 



n.l 



E 



(12) 



Using (|12|) in ( |ll| ) we are left with a new polynomial 
whose coefficients are sums of a n i x b m k and have 
to be all put to zero. With the coefficient equations in 
hand, it is easy to prove that assuming all the b m k co- 
efficient up to m + k = N are known and writing the 
equations Vi + j = (N + 1) — 3, c L j = 0, we can compute 
out of all the a n i all but three b m k with m + k = N + 1 . 
This is somewhat similar to mass reconstruction prob- 
lems from galaxy surveys where one cannot avoid the 
mass sheet degeneracy. The situation is however worse 
in our case since not only constant convergence but also 
translations and a whole class of a n k realization depen- 
dent complex deformations are indiscernible. Thus, with 
the only knowledge of the B component of the polariza- 
tion one cannot, with the first order eq. (0), recover the 
projected mass distribution. 



r coefficient 


Zgal = 1 


= 2 


EdS, Linear 


0.42 


0.60 


n = 0.3, A = 0.7, Linear 


0.31 


0.50 


n = 0.3, A = 0.7, Non Linear 


0.40 


0.59 



TABLE I. values of r, the cross-correlation between two 
source planes (z ga i and 2 cmb = 1100 ) for different models. The 
adopted filter scale (see Sect. IIIC for details) is 2 arcmin for 



both weak lensing survey and Cosmic Microwave Background 
observations. Non-linear P(k) has been computed using Pea- 
cock and Doddsmethod 



III. CROSS-CORRELATING CMB MAPS AND 
WEAK LENSING SURVEYS 

A. Motivations 

Even with the most precise experiments it is clear that 
clean detections of B component will be difficult to ob- 
tain. The magnetic polarization amplitude induced with 
such a mechanism is expected to be one order of magni- 
tude below the electric one ||isj| . Besides even if we know 
that there is a window in angular scale where the other 
secondary effects will not interfere too much with the de- 
tection of the lens-induced B |f23[ , few is known about 
removing the foregrounds |2^| to obtain clean maps re- 
construction algorithms would require. 

These considerations lead us to look for complemen- 
tary data sets to compare B with. Although the source 
plane for weak lensing surveys [|| is much closer than 
for the lensed CMB fluctuations, we expect to have a 
significant overlapping region in the two redshift lens 
distributions, so that weak lensing surveys can map a 
fair fraction of the line-of-sight CMB lenses. Conse- 
quently, weak lensing surveys can potentially provide us 
with shear maps correlated with B, but which have dif- 
ferent geometrical degeneracy, noise sources and system- 
atics than the polarization field. 

The correlation strength between the lensing effects at 
two different redshifts can be evaluated. We define r as 
the cross-correlation coefficient between two lens planes: 



(k K gal / 



(13) 



In a broad range of realistic cases (see tab. Q), r ~ 40%. 
To take advantage of this large overlapping we will con- 
sider quantity that cross correlates the CMB B field and 
galaxy surveys. Moreover, cross-correlation observations 
are expected to be insensitive to noises in weak lensing 
surveys and in CMB polarization maps. This idea has 
already been explored for temperature maps Jll|. We 
extend this study here taking advantage of the specific 
geometrical dependences uncovered in the previous sec- 
tion. 
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B. Definition of 6a and &v- 

The magnetic component of the polarization in eq. 
(0) appears to be built from a pure CMB part, which 
comes from the primordial polarization, and a gravita- 
tional lensing part. It is natural to define b, in such a way 
that mimics the AB fonction dependance, by replacing 
the CMB shear field by the galaxy one. 



(14) 



In the following, we will label local lensing quantities, 
such as what one can obtain from lensing reconstruction 
on galaxy surveys, with a g »i index. This new quantity 
can be viewed as a guess for the CMB polarization B 
component if lensing was turned on only in a redshift 
range matching the depth of galaxy surveys. The cor- 
relation coefficient of this guess with the true AB field, 

that is (^ABbj, is expected to be quadratic both in P 

and in 7 and to be proportional to the cross-coefficient 
r. 

For convenience, and in order to keep the objects we 
manipulate as simple as possible, we will not exactly im- 
plement this scheme, as it will lead to uneven angular 
derivative degrees in the two terms of resulting equa- 
tions. We can, instead, decompose the effect in the A 
and V-part. These two are not correlated, since their 
components do not share the same degrees of angular 
derivation^]. Hence, we can play the proposed game, con- 
sidering the two terms of eq. (^) as if they were two 
different fields, creating two guess-quantities that should 
correlate independently with the observed B field. Fol- 
lowing this idea we build 6a as, 



= ^7* AF* + 0(k 2 ) 



(15) 



which corresponds to the A-term in eq. (|?|). The am- 
plitude of the cross-correlation between AB and 6a can 
easily be estimated. At leading order, we have 



ABb A ) = -2eye H < 7 *7^) (AP l AP^) . 



(16) 



The corresponding V correlation is 

'AB 6 V ) = -2e«eiH (d m l k d nl ^) (d m P l d n P j ) (17) 
where we have defined 

b v = eijdk^dkPi . (18) 



*generically, a random field and its derivative at the same 
point are not correlated. 



Fig. |j shows numerical simulations presenting maps 
of first order AB, its A and V contributions and the 
corresponding guess maps one can build with a low z 
shear map. The similarities between the top maps and 
the bottom maps are not striking. Yet, under close ex- 
amination one can recognize individual patterns shared 
between the maps. This is confirmed by the computa- 
tion of the correlation coefficient between the maps, that 
shows significant overlapping, between 50% and 15%, de- 
pending correlation and filtering strategy. The calcula- 
tions hereafter will evaluate the theoretical correlation 
structure between maps given in figs. [§-b and @-g & h. 

For galaxy surveys, the amplification matrix is [fill, 



(19) 



4s»i (") - W = - / dx w g ai(x) 
Jo 

(2tt)§ 

l + cos(20 fe J sin(20 fc J 
sin(2^ fc _ L ) I - cos(2(j) k± ) 



where S(k) is the Fourier transform of the density con- 
trast at redshift z(\), w is the lens efficiency function, D 
is the angular distance, and is the position angle of 
the transverse wave- vector fcj_in the k± — (k x , k y ) plane. 
Assuming a Dirac source distribution the efficiency func- 
tion is given by 



3 DzDz^z 
w E A z ) = — ^ — 

Z Q, D z . 



(20) 



Note that the Fourier components 8{k) include the den- 
sity time evolution. They are thus proportional to the 
growth factor in the linear theory. The time evolution 
of these components is much more complicated in the 
nonlinear regime (see 
Then, \ is 

\ (a) = V(x, I, k)E(l) 5(k) ^ Kcr (l, (21) 



with the integration element defined as, 



d 3 fc d 2 l 
(2tt) 3 / 2 2^ ' 



,i[k rX +(k i D( x )+i)-5] 



(it actually depends on the position of the source plane 
through the efficiency function w(z)) and where \ stands 
for either A or V. The geometrical kernel Q Kcl ' is given 
by (using eq. (pi)) 



g^ v {T,k) = ik cos {<t> k - 



6;)sin2(0 fc 



(22) 
(23) 



This function contains all the geometrical structures of 
the A and V terms. We can write the same kind of 
equation for AB. Then, the cross-correlation is 
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(e) (f) (g) (h) 

FIG. 4. The effect of the two terms of the perturbation formula. Top row, the lens effect is the sum of the lenses up to 
recombination. Bottom row, we use the same line -of-sight mass fluctuations but only up to redshift unity, it represent our 
'local' lensing survey. The convergence fields (left panels) have been computed by slicing the z-axis and summing up the lensing 
effect in each slice. Lens-lens coupling (including departure from Born approximation) terms have been neglected, which is 
consistent with our first order approximation. The convergence in each slice has been created by using second order Lagrangian 
dynamics. The middle-left panels show the leading order contribution, the middle right the A contribution and the right the 
V one. In this example, the correlation coefficient between the two convergence maps, r is equal to 0.48 at 1.8'. The cross 
correlation coefficient between the guess map (f) and the real one (b) is 0.47. It is 0.37 between the real (b) and A (g) maps 
and goes down to 0.16 for the real (b) and V (h). 
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Xgal 



(AB\{a))=-2l Vix^l^kJ) (24) 



x ^ Kcr (f cmb ,fc cmb ) (5(k^)5(h mb )) (E(C)E(l cmh ) 



The completion of this calculation requires the use of the 
small angle approximation, 

' x 6(k ml )5(K mb )) =P(k)5 3 (k^ + K mb ) (25) 
~ P{k±) S 2 (k^ ± + k cmb± )S(k^ r + fc cmbr ) 
which implies 

fc gol = -fc cmb = k (26) 

and after the radial components have been integrated 
out, 



Xgal = Xcmb = X- 



(27) 



We also define the Ce(1) as the angular power spectrum 
of the E field, 

'E(l^)E(C b )) = C E (l)S 2 (l^ - Z cmb ) (28) 

Eventually one gets, 



2 « al f d 2 kd 2 l 

AB b\f(a) ) = -2 I d XWgalWcmb / (29) 



x c E {i)P{k)gf ev (i,k 



Then, integrating on the geometrical dependencies in 

''Ker 



5 h Kcr , we have 



ABb A (d)) =-2 J d X w gol w cmb 
dfcdZ 



- I —^kl 5 C E (l)P(k) 

(AE 2 ) (««^) , (30) 



2(2tt) 2 



and 



ABb v {a) 



x , k 3 l 3 C E (l)P(k) 



2(2tt) 2 



^((V£;) 2 )(v K -v Kg 



(31) 



implying that, ignoring filtering effects, we are able to 
measure directly the correlation between lensing effect at 
z cmb and any z gal a weak lensing survey can access. Since 
AE = AE ■ (1 + O(k)) we get, for the A type quantity, 



(AE 2 



(AE 2 



1 + 0(k)Y 
(AE 2 )-(l + 0(( K 2 ))) 



(32) 



The same holds for V. We are then able to construct two 
quantities insensitive to the normalization of CMB and 
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ABb A (a) 



(ftKgal) 



(33) 



and 



ABb v (d) 



((VEf 



(VAtga,) 2 



1 (Vk ■ Vk 

2 /Vk 2 



(34) 



1 



(Vk 2 



2 ' V V( V <.)' 

We implicitly defined ry like r but with Vk instead of k 



''v(Vi) 



v/((V K ) 2 ) ((V Kgal ) 2 



(35) 



We will see in Sect. [II D that they behave very much 
alike. This result is to be compared with the formula for 
(cos(0 g )) established in where the obtained quantity 
was going like ry/ (k 2 ). These calculations however have 
neglected the filtering effects that may significantly affect 
our conclusions. These effects are investigated in next 
section. 



C. Filtering effects 

In above section we conduct our calculations assuming 
no filtering. Obviously we have to take it into account! 
We will show here that the results obtained before hold, 
in certain limits, when one adds filtering effects. 

In the following, we consider, for simplicity, top-hat 
filters only. It is expected that other window functions 
will show very similar behaviors and this simplification 
does not restrain the generality of our results. Let us call 
W(x) the top-hat filter function in Fourier space 



W(x) = 2 



Ji(aO 



(36) 



Ji is the first J-Bessel function. We will also define Wi(x) 
a general function 



Wi(x) 



Ji(z) 



(37) 
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where Ji is the i th J-Bessel function, so that W = W\. 
Then, if X(a) is the value of any quantity X at position 
a on the sky, its top-hat filtered value can be computed 

as 



^ X k W(k9) e*' 5 , (38) 

Z7T 



where X is X Fourier transform. In the following we will 
note X(a} the filtered quantity at scale 9. 



The tricky thing for (ABb^J is that the CMB part and 

the low-redshift weak lensing part are a priori filtered at 
different scale. For AB, which is a measured value, its 
pure CMB part and its weak lensing part are filtered at 
the same scale 9. Hence, B reads, 



AB(a) 



V( x J,k)E(l)6(k) 



Ql cr (lk A 



->Kcr 
'A 



I, fej 



W (\k ± D + l\e) (39) 



A contrario \ is a composite value. The CMB part is 
still filtered at 9 whereas the weak lensing part (which 
comes from a weak lensing survey of galaxies) is filtered 
independently at another scale which we denote 9 gai . It 
implies that, 



(40) 



\[a) {e) =-2j V( X ,l,k)E(l)5(k) 
x£f er (zX) W(kD6 gBl )W(W). 



Taking filtering into account, the cross-correlation co- 
efficient becomes, 



< AB/(j\ h 



(0) °ti(0,Vi) 



= -2 



d% w go ,w cn 



(41) 



^0c E (i)p(k)g^(i,k) 
xw(kB e^w{w)w{\m + i\e). 

It can be shown (from the summation theorems of the 
Bessel functions) that, 



Wx(\kD + j\6) = 

-J2iW i (kBe)W i (W)(-l) i ^ 



(42) 



sini 







It is then possible to break the W ^JfcD + Zj#J into a sum 

of Wi(kD9)Wi(W) with coefficients that depend on the 
geometrical properties of our problem. Integrating over 
the geometrical dependencies of G\T eT , leaves us with only 
a few non vanishing terms in our sum, 



d<t> sin 2 (20) 



smuc 



sm < 

for the A-term and 



7r i = 1 or i 
elsewhere 



(43) 



d(f) cos</>sin 2 (20)sin(i0)sin</> 



tt/2 1=1 

3tt/4 i = 3 

7r/4 i = 5 

elsewhere 



(44) 

for the V-term. Each term can be computed exactly, 
and it turns out that the terms built from Wi, i > 1 
are always negligible compared to the ones coming from 
W\. It implies that we can safely ignore the W 3 and W 5 
in both A and V expressions, therefore it is reasonable 

to assume that w(\kD + l\ej = W{kB9)W(W). It is 

expected that other windows, in particular the Gaussian 
window function, share similar properties. Then, taking 
into accounts the filtering effects, the equations for the 
cross-correlations reduce to 



A-B(e)&A(e,6> gal ) 



AEf e) ) {K (e) K M o^ } ) (45) 



and 



A£?( e )6v( 



V£f fl) )(VK (e) VA 



so that our correlation coefficients can be written, 



X 



A(e,Vi) 



~ r (0,e g ai), 



\0) 



(46) 



(47) 



and 



v(e,e g ,i) - -2 rv ( e ' e sai) 



v 4)> 



(48) 



The results obtained in eqs. d33fj34j ) are thus still formally 
valid. Actually, eqs. ( f47jj48| ) simply tell that filtering 
effects can simply be assumed to act independently on 
the lensing effects and on the primary Cosmic Microwave 
Background maps. We are left with two quantities that 
only reflects the line-of-sight overlapping effects of lensing 
distortions. 



D. Sensitivity to the cosmic parameters 

We quickly explore here the behavior of X^ in differ- 
ent sets of cosmological parameters. These quantities 
only depend on weak lensing quantities. Ignoring the 
VLq dependence in the angular distances and growing fac- 
tor, one would expect (k 2 ) to scale like Oq. Yet, be- 
cause of the growth factor, the convergence field exhibits 
a weaker sensitivity to O . Assuming A = and a power 
law spectrum, we know from |16| that («; 2 al ) oc fig 66 
for z gal = 1. The same calculation leads to (K cmb K gal ) oc 
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E. Cosmic variance 




0.2 0.4 0.6 
Q n value 



1.0 



FIG. 5. («(e)« AB i(e Ral )) / \ K g a i(9 gal )/ f° r a CDM model con- 
sistent with the values of (fio, A). 8 — 8 &al — 2'. 




0.2 0.4 0.6 
Q n value 



1.0 



FIG. 6. (V« W • V Kgol(9gal) )/((VK gal(£)gal) ) 2 ) for a CDM 
model. 9 = gal = 2'. 



)1.915 



ill 68 , ((VK gal ) 2 ) oc nl m and (v«; cmb • Vk 8s1 | oc CI 

Then, in this limit, the quantities X\ have a very low 
dependence on flo : 



0.005 



X A oc and X v cx ^ 

Eventually, the X^ quantities should exhibit a seizable 
sensitivity to A; changing A increases or reduces the size 
of the optic bench and accordingly the overlapping be- 
tween K cmb and K gal . 

Figs. H and || present contour plots of the amplitude of 
X A and X v in the (fi 0l A) plane for CDM models. They 
show the predicted low Qq sensitivity and the expected 
A dependency. Both figures are very alike. This is due to 
the fact that the dominant features are contained in the 
efficiency function dependences on the angular distances. 



In previous sections we looked at the sensitivity of ob- 
servable quantities which mixed galaxy weak lensing sur- 
veys and CMB polarization detection. It is very unlikely 
that both surveys will be able to cover, with a good reso- 
lution and low foreground contamination, a large fraction 
of the sky. It seems however reasonable to expect to have 
at our disposal patches of at least a few hundreds square 
degrees. The issue we address in this section is to es- 
timate the cosmic variance of such a detection in joint 
surveys in about 100 square degrees. 

The computation of cosmic variance is a classical prob- 
lem in cosmological observation J2(J . A natural estimate 
for an ensemble average (X) is its geometrical average. 
If the survey has size E then, 



X= i / d 2 aX(a) 



(49) 



For a compact survey with circular shape of radius H we 
formally have, 



X= f ^X(k)W(kZ). 



(50) 



For sake of simplicity this is what we use in the follow- 
ing but we will see that the shape of the survey has no 
significant consequences. 

Taking X as an estimate of (X) (the ensem- 
ble average of X), leads to an error of the order 

i^ 2 ^ — (X} 2 which usually scales like if the sur- 

vey is large enough. 

When we are measuring X\ on a small patch of the 
sky, we are apart from the statistical value by the same 

kind of error. We can neglect the errors on lAE' 



(\7E) 2 J, ((VK gal ) 2 )and («: gal ); those may not be the 

dominant source of discrepancy and can even be mea- 
sured on wider and independent samples. The biggest 

source of error is the measure of (ABb^ \. It is given by, 



a = 



ABbh - ABbi 



— (ABbh — ABh 



(51) 



The computation of ( |5l| ) is made easier if we write ex- 
plicitly the geometrical average as a summation over N 
measurement points (TV can be as large as we want), 



X 



1 N 



(52) 



we then developed (pT[) , and replace the ensemble av- 
erage of the summation sign by the geometrical average 
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over the survey size. We are left with a sum of correlators 
containing 8 fields taken at 2, 3 and 4 different points. 
The calculations can be carried out analytically if we as- 
sume that all our fields follow Gaussian statistics, which 
is reasonable at the scale we are working on. In that case 
indeed, we can take advantage of the Wick theorem to 
contract each of the 8 fields correlators in products of 2 
points correlation functions. By definition, ( |5l] ) contains 
only connected correlators, moreover the ensemble aver- 
ages (AB) and (fy) vanish, therefore only a small fraction 
of correlators among all the possible combination of the 
8 fields survive. We can use a simple diagrammatic rep- 
resentations to describe their geometrical shape. All the 
non vanishing terms in C§ are given in Fig. [jj Each 
line between two vertex represents a 2 points correlation 
function such as (X(ai)X(a2)) , and the different sym- 
bols at the vertex correspond to different X fields (the 
cross stands for AP, the dot for 7 cmb , and the open dot 
stands for 7 ga ,). The .A-terms represent terms where the 
two top (and the two bottom) AB and are taken at 
the same point, but top and bottom fields are not at the 
same place. The i3-terms are three points diagrams: the 
top AB and are at the same point whereas the right 
and left bottom vertexes are at two different locations. 
The C terms are four-points diagrams, where each vertex 
is at a different point. To illustrate our notations, let us 
write B\ c as an example, 

#L = (7cmb(<3i)7 gal (a 2 ))(7 g! .i(a3)7cmb(<3i))x 
(AP{d 1 )AP(a 1 ))(AP(d 2 )AP{d 3 )) 

We only focus on the calculation of the A terms be- 
cause we can use the approximation that 



2 x 




+ ■ ! + A 

II / \ 

:0 JLQ 


W 3 jfi 



+2 x 



? 

! +2x 



V 

/ \ 



A 



3 c 



-2 x 



\ i \ v ' 



o 



o • 



+2x1; \ )+2 x ;<\ 



„4>£>C. 



(53) 



Indeed, in perturbative theory, if the survey is large 
enough, the n-points correlation functions naturally dom- 
inates over the n + 1-points correlation function. This is 
true as long as the local variance is much bigger than 
the autocorrelation at survey scale and we assume the 
surveys are still large enough to be in this case. 
The general expression for any A diagram is 



A 



-2 x 



Z>(X«nbl, Ubl, *«nbl) ^(Xc mb 2, Z cmb 2, K mh2 ) (54) 
^(Xgall, Zg.ll, fcgall) f(Xg.l2) ig-l2l fcgal2) 



FIG. 7. Diagrammatic representation of the terms con- 
tributing to the cosmic variance of the correlation coefficients. 
In this representation the vertex x • represents AB; the cross 
stands for the AP part, the dot for 7 cmb . The other vertex 
xo represents any \\ the open dot stands for 7 g „i. The solid 
lines connect AP terms and the dashed ones the 7-s 



x^ Ker (/c„bi, fc cm „i J ^ Kor (Z cm b2, K mb2 J 

Xgf er (C,l,fe g a l lJ^ Ker (C2,C2x)M\- , 

xW(\K mhl± B + f mbl |0) W(\K mh2± V + l ctnh2 \6) 
xW{K M J>iO sal )W{l^x6) W(k^ 2± D 2 9^)W(l^ l2 6) 

xW (|fc g »iuDi + Ci + LuDi + C*i|s) 
xW (|*Uu D 2 + f gal2 + k cmh2± D 2 + Cabals') 
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where Mi gives the 2-point correlations associated with 
the lines of the diagram. For example : 



(55) 



We explicit in the following the computation of A\. 
The other terms follow the same treatment or can be ne- 
glected. The lines in the A\ diagram give us the relations 



Mi ={5(k^ a )6(h mbl )) ^(fc ga , 2 )<5(fc cmb2 
x (E(l sall )E(l S!d2 )) (E(l cmhl )E(l cmh2 ) 



^cmbl ^gall — fcl 

fccmb2 = -fcgal2 = k 2 

Lmbl = ~lmb2 = Lb 

^gall = ~^gal2 = ^gal 



(56) 



Then, using these relations and the small angular ap- 
proximation, we have : 

A\ = 4: dxidx2 w cmb iw gal iWcmb2W gal2 (57) 

d 2 kxd 2 k 2 d 2 Z gal d 2 / cmb 
~(2^Y i (2tt) 4 
xC E (l^)C E (l atab )P(k 1 )P(k 2 ) 
*9f c *(i mh M)Qf cr (^k 2 ) 

x^ Kcr (C„-fci)^ Kcr (-r g a„-fc2) 

xWflfaD + T mb \0)W{\k 2 T> + t mb \6) 
xW{k{D6 s S)W{k 2 m sal )W 2 {l &ai 6) 

We apply the decomposition of W\ (^\kT){x) -\-f\0\ we 

used in eq. (^2|). The geometry of our problem is 
the same and the result (Q) still holds for the terms 

in Wi(|fciD( X i) + f cmb |0) and W t (|£ 2 D( X2 ) + C b \6) , 

This however is not true for (j? ga i + L„,h\^j for which 
the application of the re-summation theorem does not 
bring any simplification. Then, neglecting all the W3 
parts and after integration on the cfiki, for the A-term, 
we have, 

/gal 
dXldX2 Wc„blW g ,,lW cmb 2W g al2 

dfcidfc 2 d 2 Z ga ,d 2 Z cmb 4 4 

x C B (Z gal )C B (/ cmb )P(fci )P(fe) (58) 
x VF 2 (|f ga , + f cmb |S) cos 2 2 (0, cmb - ignl ) 
x W (fciD0 gal ) VF(fc 2 D# gal ) (Z ga! 0) 
xW^(fciD6»)Vy(fc 2 D6')VK 2 (Z cmb 6i). 

Note that for the evaluation of the V part, using the same 
kind of method, we obtain the same equation as eq. ( |5^ ) 
where Z 4 al ^ mb is replaced by I 2 J 2 mh k\k 2 / 2 . 




Filtering scale (arcmin) 



FIG. 8. Comparison between y/2 Af /signal A (solid line) 
and yj2 Af /signal v (dashed line). The C e are from a 
Q. = 0.3, A = 0.7 model. The survey size is 100 deg 2 , and 
Gaussian filters were used. 



We can get rid of the remaining W 2 Hi gal + / cmb |^ 
with another approximation. The power spectrum C E (l) 
favors large values of I whereas W (|/ gaI + Z cmb |S) will be 
non-zero for |Z gal + Z cmb | ~ 1/H. Then for typical survey 
size of about one hundred square-degrees, |2 go i+/ omb | <C 1% 



and we can assume Z gal 
limit, cos 2 2 Ui, - <j>i 



-Z cmb and Z g ai + Zc 



e. In this 



< = 



Idl 

1 

dX^cmb 



l 8 C 2 E (l)W 4 {W) 
kdk 



1 and A\ can be written 
d 2 e 



2tt 



rW({e 5) (59) 



2tt 



p(k)w{km)w{km^) 



which is essentially the cosmic variance of (AE 2 ), for the 
A part and of ((V.E) 2 ) for the V one (where I 8 in eq. 
is replaced by l 4 k 2 k 2 /2). Finally we have, 



At _ 2^ fdll 9 C 2 E (l)W?(lO) 

: %)&A( Mg a,)) 2 _ ^{fdll5C E (l)W}(lO)f 

oc Cosmic variance of Ai? 2 



(60) 



where £ = 7rS 2 in case of a disc shape survey. We show in 
Fig. H numerical results for a 100 deg 2 survey although 
the numerical calculations were done with a Gaussian 
window function instead of a top-hat. 
Numerically, for 8 — 10', we get 



A 



(3.7%) 2 

<b w w,)) 2 ~ s / 100 deg2 



(61) 



We expect that for the same reasons, the A 2 terms will 
be dominated by the weak lensing variance. Yet a correct 
evaluation here is harder to reach. We have made this 
estimation within the framework of a power law P{k). 
With this simplification in hand, we can write for A 2n 
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= 5', 6» gol = 2.5' 
= 5', gol = 5' 
> = 10', 6> aal = 5' 



CosVar ((k 2 )) CosVar (((Vk) 2 )) 

Q = 0.3 n = 1 = 0.3 Q = 1 

2.94% 1.86% 2.88% 2.07% 

3.02% 1.87% 2.23% 1.75% 

3.54% 2.03% 4.25% 3.02% 



TABLE II. Values of the cosmic variance of < K 2 ) and 
^(Vac) 2 ^) for different models and different filtering radius. 
The size of the survey is 100 deg 2 . For the ft = 0.3 (fio = 1) 
model, we use 5 (7) independent ray-tracing realizations (see 
[pif) to estimate the cosmic variance in a 9 deg 2 survey, which 
is then rescaled to the cosmic variance we should obtain for 
a 100 deg 2 survey. Given the low number of realizations, the 
values here can only be used as a good estimation of the or- 
der of mag nitude of CosVar ((k 2 )) and CosVar (((Vk) 2 )) . 
It also seems, from these figures that the cosmic variance of 
{^(Vfi;) 2 ^ is more degraded by the difference in filtering beams 
than the other. 



(we focus only the A part, but the same discussion holds 
for the V observable.) 



A A 



(-B(0)&A(0,0 gal )) 

' /, d 2 fc 1 d 2 fc 2 P(fc 1 )P(fc 2 )cos 2 (0 fcl -<f> k2 ) 
< Wj»(M)Wf»(M,.,)W?(|fci + k 2 \E) 



(62) 



(63) 



The last integral behaves essentially like the cosmic 
variance of (k 2 ). More exactly, it goes like I/a/2 this 
variance. It should even be smaller, because of the ex- 
tra cos 2 factor. We evaluated this cosmic variance using 
the ray-tracing simulations described in |2^| . These sim- 
ulations provide us with realistic convergence maps (for 
the cosmological models we are interested in) with a res- 
olution of 0.1', and a survey size of 9 square degrees. 
The sources have been put at a redshift unity, and the 
ray-lights are propagated through a simulated Universe 
whose the density field has been evolved from an initial 
CDM power spectrum. The measured cosmic variance of 
(k(6)K(9 is about 3% (see Table ||) when filtered at 
scales 6> gal = 5' and 9 = 10' for a flo = 0.3 cosmology. An 
estimation of A A n is then given by, 



A A 



2.12% \ 2 1 
r J E/100 deg 2 ' 



Since ry is very comparable to r, we very roughly esti- 
mate Aj n 



A v 



1 



2.12% y 

r J E/100 deg 2 ' 



(65) 



The same considerations gives 



= 5', 6» gal = 2.5' 
= 5', gal = 5' 
' = 10', 6» , = 5' 



CosVar (X A ) 

n = 0.3 fi Q = i 

6.44% 4.77% 
6.58% 4.79% 
8.71% 6.73% 



CosVar (X v ) 
Qo = 0.3 fio = 1 
6.06% 4.72% 
4.99% 4.23% 
9.49% 7.62% 



TABLE III. Values of the cosmic variance of X§. The sur- 
vey size is 100 deg 2 . We used the results presented in Table^] 
and Fig. @. The r\ parameters are assumed to be equal and 
set to 0.4. We didn't take into account the filtering effects in 
the definition of r. The difference due to filtering correction is 
small, though. From this estimations, we can expect a cosmic 
variance for X^oi less than 10% for realistic scenarii. 



A 



(2.12%) 2 



(66) 



There is no r dependency here; the diagram cross- 
correlates K cmb and K gal . 

We can approximate the remaining .A-terms. They 
should be smaller than the former. We have 



A 



h 1 (2.12%x3.7%) 2 /D x2 
3n ~ 12 r/1m , 2 — \ B {8) b Wfi^)J 



ri E/100 deg 



«-4 



and 



(2.12% x 3.7%) 2 
E/100 deg 2 



Then, only the A\ and A% terms (boxed on Fig. 0) con- 
tribute substantially to the cosmic variance of X\ . Since 
A\ and A% are respectively the cosmic variance of (A.E 2 ) 
(resp. ^(VB) 2 ^) and of (k 2 ) (resp. ^(Vk) 2 J>), we can 
write the variance of Xu as 



CosVar(A- A ) = 
CosVar ((AE 2 )) 



(67) 



2r 2 



CosVar ((k 2 )) 



and 



CosVar(A- v ) - 
CosVar (((VE) 



2r% 



^ ) CosVar 



Table presents numerical results for various filtering 
scenarii and models. 

The two quantities, 6 a and fry, lead to similar cosmic 
variance that are rather small. Obviously it would be 
even better to use b = &a + 6 v • For such a quantity the 
resulting cosmic variance for the cross-correlation coeffi- 
cient should even be smaller, by a factor v2, although 
a detailed analysis is made complicated because of the 
complex correlation patterns it contains. 
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IV. CONCLUSION 

We have computed a first order mapping that de- 
scribes, in real space, the weak lensing effects on the 
CMB polarization. In particular we derived the explicit 
mathematical relation between the primary CMB po- 
larization and the shear field at leading order in lens 
effect. It demonstrates that a i?-component of the po- 
larization field can be induced by lens couplings. We 
have shown however that the -B-map alone cannot lead 
to a non-ambiguous reconstruction of the projected mass 
map. 

Nonetheless, the J5-component can potentially exhibit 
a significant correlation signal with local weak lensing 
surveys. This opens a new window for detecting lens 
effects on CMB maps. In particular, and contrary to 
previous studies involving the temperature maps alone, 
we found that such a correlation can be measured with a 
rather high signal to noise ratio even in surveys of rather 
modest size and resolution. Anticipating data sets that 
should be available in the near future, (100 deg 2 survey, 
with 5' resolution for galaxy survey and 10' Gaussian 
beam size for CMB polarization detection) , we have ob- 
tained a cosmic variance around 8%. Needless is to say 
that this estimation does not take into account system- 
atics and possible foreground contaminations. It shows 
anyway that Cosmic Microwave Background polariza- 
tion contains a precious window for studying the large 
scale mass distribution and consequently putting new 
constraints on the cosmological parameters. 

In this paper we have investigated specific quantities 
that would accessible to observations. They both would 
permit to put constraint on the cosmological constant. 
The simulated maps we presented here are only of il- 
lustrative interest. We plan to complement this study 
with extensive numerical experiments to validate our re- 
sults (in particular on the cosmic variance), and explore 
the effect of realistic ingredients we did not include in 
our simple analytical framework, a shear non-gaussianity, 
lens-lens coupling and so forth. 
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